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Abstract. We modify the scalar Einstein-aether theory by breaking the Lorentz invariance 
of a gravitational theory coupled to a Galileon type scalar field. This is done by introducing a 
Lagrange multiplier term into the action, thus ensuring that the gradient of the scalar field is 
time-like, with unit norm. The theory can also be considered as an extension to the mimetic 
dark matter theory, by adding some derivative self interactions to the action, which keeps the 
equation of motion at most second order in time derivatives. The cosmological implications 
of the model are discussed in detail. In particular, for pressure-less baryonic matter, we show 
that the universe experiences a late time acceleration. The cosmological implications of a 
special coupling between the scalar field and the trace of the energy-momentum tensor are 
also explored. 
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1 Introduction 

Lorentz invariance seems to be a fundamental property of nature which has been tested 
within the context of the standard model of particle physics with a very high accuracy [1]. 
However, in the gravity sector, Lorentz invariance cannot be tested to such a high precision. 
All the gravitational experiments put a mild bound on the validity of Lorentz invariance 
[2]. In the dark sector of the universe, we have some constraints on the breaking of Lorentz 
invariance in dark matter [3] and dark energy sectors [4]. For a very good review, see [5]. This 
means that one can build a Lorentz invariant theory which breaks the Lorentz invariance but 
satisfyie the experimental data. In this sense, it would be plausible to have a gravitational 
theory, describing the large scale structure and dynamics of the universe which breaks Lorentz 
invariance. 

Another motivation for considering theories that break Lorentz invariance comes from 
quantum gravity. A very interesting example of such a theory is the Horava-Lifshitz gravity 
where adding higher order Lorentz violating terms to the action together with different scaling 
dimensions for space and time makes the theory power counting renormalizable at the ultra¬ 
violet level [6]. In the infrared region, in the Horava-Lifshitz theory, the scaling dimension 
of time becomes unity. Therefore, the Horava-Lifshitz theory can in principle be written in 
a generally covariant form [7]. However, due to a fixed foliation of space-time into space-like 
hypersurfaces, the theory should break Lorentz invariance dynamically. 

One of the best candidates for a generally covariant gravitational theory with a preferred 
time-like direction is the Einstein-aether theory, proposed in [8]. The preferred direction can 
be imposed to the theory by introducing a time-like vector field to the action through a 
Lagrange multiplier. The action of the theory then becomes 
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( 1 . 1 ) 






( 1 . 2 ) 


where A is a Lagrange multiplier and the tensor is defined as 

K Z = cog^gxa + ci« + C 2 5^5" + c 3 u^u u g X(T , 

where Cj, i = 0,1, 2,3 are the dimensionless free parameters [ 8 ]. The stability of the theory has 
also been studied and it turns out that the theory has a non-empty region in the parameter 
space in which no instabilities occur [9]. The PPN parameters were also obtained [10] and 
it was shown that all PPN parameters are identical to the standard GR if one imposes the 
conditions c 2 = (—2 cf — c\c 3 + c|)/3ci and C 4 = —c^/ci respectively, for a discussion of this 
issue see also [ 11 ]. 

The other interesting fact about the Einstein-aether theory is that it can be considered 
as a covariant version of the Horava-Lifshitz theory in the infrared sector. In fact, one can 
prove that the above theory is identical to the non-projectable Horava-Lifshitz theory for 
a special case of the aether field, u^ = NX7^T, where N is a normalization factor [12], 
However, the action (1.1) can only produce the lowest terms of the Horava-Lifshitz theory 
and the search for a generally covariant theory, which produces higher order terms of the 
Horava-Lifshitz theory is still a subject of research. 

It will also be interesting if one could write an action for a scalar field with a preferred 
time-like direction. In this case the covariant derivative of the scalar field can play the 
role of a time-like vector. In [13] the authors have suggested such a gravitational model 
by substituting u^ by in the action (1.1). However, this choice separates only the 
longitudinal mode of the vector field and also implies that has a unit norm. It was 
shown in [14] that the theory proposed in [13] is identical to the projectable Horava-Lifshitz 
gravity which has instabilities [15, 16]. The problem lies in the aether action itself, where 
substituting the aether field with V M <(> produces some fourth order derivative terms. These 
higher order derivative terms will produce ghost instabilities to the action. Hence, in order to 
write a scalar action which breaks the Lorentz invariance dynamically, one should introduce 
a kinetic term for the scalar field which has no instabilities. 

One of the main motivations of the present paper is to write down an action for a 
“healthy” scalar-aether theory. We first note that in writing the action (1.1), the authors 
have introduced the most general canonical kinetic term for a vector field. In the scalar case, 
the canonical kinetic term is of the form d^d^cj), which is the same as the term used to break 
the Lorentz invariance. In fact, one adds a term X^d^d^cj) + 1) to the action where A is a 
Lagrange multiplier. In this paper, we want to add some higher order kinetic interactions to 
the theory in a way that the resulting theory has no Ostrogradski instability. In this sense 
one should add some healthy higher order derivative interactions to the action. Such terms 
are known as Galileons [17]. 

The Galileons are scalar fields where having even higher than second order time deriva¬ 
tives in the action produce second order field equations. The Galileons were originally written 
in flat Minkowski space [17] where they have an additional symmetry 

<f> ->• 4> + + c, (1.3) 

which is lost in the covariant version [18] 1 . In order to have a covariant a Galileon theory, one 
can replace the partial derivatives with the covariant derivatives. However, because of the 

x The most general theory of two graviton degrees of freedom (dof) plus a scalar dof writing in terms of 
curvature invariants and the covariant derivatives of the scalar field was first introduced by Horndeski in [19], 
and independently rediscovered in [20]. The covariant Galileon theory which we are using in this paper is 
then a special case of Horndeski theory. 
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non-commutativity of the covariant derivatives, one obtains some higher order terms which 
produce Ostrogradski instabilities in the curved background. In order to remove these non- 
healthy terms one should add to the action some new terms which do not respect the Galilean 
symmetry (1.3). This will produce a covariant Galileon action [18]. One can also generalize 
the Galileons by substituting the coefficients of each term with an arbitrary function of the 
scalar field [ 22 ]. 

Recently, a very interesting model for dark matter sector of the universe was proposed in 
the literature, known as the “mimetic dark matter” [23]. The model has originally suggested 
a new action for a conformally invariant gravitational theory. In this sense, the authors have 
added a scalar field to the theory and built the gravitational action using the effective metric 

9/jtv = {g^d^d^g^. (1.4) 

The theory is obviously invariant under the change g IJV —>• VL 2 g^ v , with Q an arbitrary func¬ 
tion. One can easily observe that in this theory the scalar field cj) is subject to the constraint 
g^ u d^(j)d v (j) = — 1, which implies that the vector u^ is tinre-like [23, 24], So the vector 

field u^ will be tangent to time-like geodesics, and will correspond to the 4-velocity of a dust 
particle. This geometrical dust corresponds to dark matter. The mimetic theory can then 
explain the dark matter content of the universe as a consequence of the constraint equation 
g^d^duc/) = — 1. This also happens in the scalar Einstein-aether theory where one adds a 
similar constraint to the action through a Lagrange multiplier. This suggests that the model 
introduced in the present paper can be considered as an extension of the mimetic dark matter 
theory with the addition of the most general higher derivative self interaction which limits 
the equations of motion to second order. Such a modification makes dark matter sector of 
the universe imperfect [25]. We should note that one can generalize the constraint equation 
by rewriting it as g^d^d^cf) = —fi(cj)) 2 [26]. This condition implies that the vector field 
U[i = dfj,(j)/g is tangent to time-like geodesics, corresponding to the 4-velocity of dark matter. 

The present paper is organized as follows. In Section 2 we introduce the action of the 
Lorentz violating Galileon theory and obtain the equations of motion. In Section 3 we study 
the cosmological implications of the theory, showing that it can explain the accelerated 
expansion of a matter dominated universe, having only zero pressure dust in the energy- 
momentum of the universe. In Section 4 we add a special interaction between the aether field 
and the trace of the energy momentum tensor of ordinary matter and study the cosmology 
of this modified gravity model. Last Section is devoted a discussion of our results. 


2 Generalized scalar Einstein-aether theory 


We propose a generalization of the scalar Einstein-aether theory based on a gravitational 
action of the form 


S 



K 2 R + 03 X 3 + 04 X 4 + 05/25 + + 1) + C m , 


( 2 . 1 ) 


where C m is the matter Lagrangian and <= V M </>. In the above action we have dropped 
the tadpole Galileon term and also absorbed the quadratic term into the term A ((j>u,<i>^ + 1). 
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We have also introduced the terms £j, i = 1,2,3, defined as [18] 
£3 = 

£4 — (4 > a4 


1 


£3 — 


(□0) 3 - + 2cj) ^ " - 6^^G i 


vpV 


( 2 . 2 ) 

(2.3) 

(2.4) 


The action (2.1) is similar to the action for mimetic dark matter with derivative interaction 
terms for a scalar field [24]. 

Variation of the action with respect to A gives 


4>u^ = -1, 


(2.5) 


which, as was mentioned before, defines a preferred direction for the space-time, and dynam¬ 
ically breaks the Lorentz symmetry of the theory. After using Eq. (2.5) and its derivative 


= °> 

we obtain the reduced form of the metric field equations as 

v + a-i T^ u + aq T* v + as T^ u — A 4>p<t> v , 
where T'^u is the ordinary matter energy momentum tensor defined as 

rp _ —2 5{^/—gC m ) 
sTa ’ 

and we have defined 


( 2 . 6 ) 


(2.7) 


T nv = 

T U = — 2</vD</> + 24 >p P 4 ) p v + ((D0) 2 - 4)p<rV a ) {dpu - 2<^</>„) 


( 2 . 8 ) 


— + 2 4> p {R pll 4>i/ + Rpv4> p ) + — 2^>p7? pCT 0 (T 5 /i j / + 2cj) p (f) a Rpp Va , (2.9) 


and 


T pu = -3((D</>) 2 - ^ctA^ <TA )^ + - ~4>p4>uR) + 3(j> ,T n(t>(R' Tll (/> I/ + R a u4>p) 

+ 3D (jxjf^RpavX - 6 (f)po4> ap( !>pv + ‘Z(j)o\R aX (t>p,(t ) v - 3(/) a R CTX ((j)xp(j>v + 4>\v4>p) 

— 3 4> a 4> Xp (RpXaprf’v + RvXap^p) + 3 4> a 4> X (RpaXpfit 1 v + Rf<j\p4 }P p) + 3(j) (T (l)x4>pK.R <JpXK 9pv 

+ ((D0) 3 - 3D ^paVl + 2 4>pa^ X ^x P ^ (9pu - Mu) + 3((t>„R aX (t>x)(M - D^). 

( 2 . 10 ) 

The aether field equation of motion can be written as 

03^3 + 04^4 + as£s — 2 V m (A(/> m ) = 0 , ( 2 . 11 ) 
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where we have defined 


£ 3 = -2 ( ncj )) 2 +2 (2.12) 

£* = - 4 P -/ 0 3 - - 2 {u<i>)r 

+ 8 + 4<\> pv R pv - , (2.13) 


and 

f 5 = - 2(Dp 4 + 3(D<))) 2 ( 4 ^^" - R + 24>pR» v (t>^j - **) + 12(D«/»)^^ 

-12 (n<i>)<i> lt <i> v <i> piT RW' a - +12 + 3 ^r v )R 

- 6 (^(TK^pR^M - 12<l>v p R pa (t>a U - Hv P 4>a\R V(JpX + 12 Mvtpar xR^ 

+ 3 (c^ u R u y p )R - G(t> v R pp R pa r - WA p R aX R uapx + 3 ^ P R V anX R p(TKX . (2.14) 

Note that the energy-momentum tensors T pu have the property that their covariant 
derivatives become proportional to E % [17] 

V # ‘7i / = ^i, i = 1,2,3. (2.15) 

So, the covariant derivative of the equation (2.7) is reduced to 

V^Tp = -^4> u [a 3 £ 3 + a 4 P + a 5 £ 5 - 2V M (A<jp]. (2.16) 

The energy-momentum tensor of ordinary matter is then conserved if the equation of motion 
for the aether held is satisfied. In the following we will assume that the energy momentum 
tensor of ordinary matter is conserved X7 p T pu = 0 . 


3 Cosmological implications 

In this Section we study the cosmological implications of the Lorentz violating Galileon 
theory. We will restrict our study to homogeneous and isotropic cosmological models with 
the line element given by the flat Friedmann-Robertson-Walker metric 

ds 2 = — dt 2 + a 2 (t)dx 2 . (3.1) 

We also assume that the aether held and the Lagrange multiplier are homogeneous and 
therefore have the form (f> = 4>(t) and A = A(f), respectively. We also assume that the matter 
energy-momentum tensor, describing the matter content of the universe, has a perfect fluid 
form 


T p v = diag[ — p{t),p{t),p{t),p{t )\, 


(3.2) 


where p(t) is the matter energy density, while p(t) represents the thermodynamic pressure. 
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3.1 The cosmological evolution equations 


Due to the choice of the geometry as homogeneous and isotropic, the constraint equation 
Eq. (2.5) completely determines the aether held <f> as 

(j> = t + c i, (3.3) 

where c\ is an arbitrary integration constant. For a homogeneous and isotropic cosmological 
model, Eq. (3.3) is obvious since the FRW metric already has a time-like preferred direction 
d/dt, and the aether vector should be identical to that direction up to a shift. 

The metric and scalar held equations can then be obtained from Eqs. (2.7) and (2.11) 
as 

^(15a4 + 2 k 2 )H 2 — 21a5 H i — 3a 3 H — p + A = 0, (3-4) 

(3a 4 + 2k 2 — 605 H)H + ^(3«4 + 2 k 2 )H 2 — 6 a 3 H 3 + p = 0, (3-5) 


and 


6(15 a 5 H 2 - 12 a A H + a 3 )H + 90 a 5 H 4 - 108 o^H 3 + 18 a 3 H 2 - 2(3A H + A) = 0. (3.6) 

In Eqs. (3.4) - (3.6) H = H(t) = a{t)/a{t) denotes the Hubble parameter. 

One can solve the aether equation for the Lagrange multiplier to obtain 

A = 3H(a 3 — 604 H + 5 a 3 H 2 ) -\ —(3-7) 

a 6 

where C 2 is an arbitrary constant of integration. Substituting Eq. (3.7) into Eqs. (3.4) and 
(3.5), one obtains 

-(2k 2 + 3a 4 )H 2 — — p = 0, (3-8) 

2 a A 

(3a 4 + 2 k 2 — 6a 3 H)H + ^(3ct4 + 2 k 2 )H 2 — 6 a^H 2, + p = 0. (3-9) 

One can easily see that a 3 introduces higher order Hubble parameter to the Friedmann 
equation and aq modifies the gravitational constant of the theory. Solving Eq. (3.8) for p 
gives 


p = ^iL 2 (3«4 + 2 k 2 — 4«5 H) + (3.10) 

.z a 

One of the most common equations of state, which has been used extensively to study 
the properties of compact objects at high densities is the linear barotropic equation of state, 
with p = cop with u = constant € [0,1]. Assuming the equation of state of the ordinary 
matter as having a linear barotropic form and substituting Eq. (3.10) to Eq. (3.9), we obtain 
the basic cosmological evolution equation of our model as 

(2k 2 + 3«4 — 6a 3 H)H + ^(u; + 1) (2 k 2 + 3«4 — 4 a§H) H 2 + = 0 . (3.11) 

Z CL 
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3.2 The general dynamics of the SEA cosmological models 

In order to obtain a simpler form of Eq. (3.11) we rescale the parameters 04 , 05 and C 2 so 
that 


2kt k z 

a 4 = — m, a 5 = yn 4 , c 2 = 2k 2 s 4 , 


where m, n 4 , si are constants. Then Eq. (3.11) takes the form 


(1 + m — n\H)H + - (oj + 1) 


1 + m -ni 77 

3 


H 2 + ^ = 0. 


(3.12) 


(3.13) 


We introduce as an independent variable the redshift z, defined as 1 + z = 1/a. Therefore 

(3.14) 


dH dHdz , , TT dH 
— = ——- = -(1 + z)H——. 

dt dz dt dz 


Moreover, we represent H(z) as H(z) = Hoh(z), where Hq is the present value of the 
Hubble parameter. By rescaling again the coefficients ni and si so that 


ni = 


n 


s = 


si 


(3.15) 


we obtain the basic evolution equation of the Hubble parameter as 


(1 + z)h(z) [1 + m — nh(z)\ ^ ^ {u + 1) 


1 + m— —nh(z ) 


h 2 (z) + ws(l + z) 3 . (3.16) 


By rescaling the density as p(z) = 2n 2 HQr(z) one obtains 


r(z) = 


1 + m — —nh(z ) 


h 2 (z ) + s(l + zf. 


(3.17) 


By rescaling the constant 03 as 03 = 2k 2 uHq, where u is a constant, and the Lagrange 
multiplier A as A = 2n 2 H/jA, we obtain 


A (z) = 3 h(z) 


u-2mh(z)-\ — nh 2 (z ) +5(1 + ^)" 
6 


(3.18) 


As an indicator of the accelerated expansion we introduce the deceleration parameter, 
defined as 


d 1 i _ H(t) 
dtWfj ~ 1 ~ “ H 2 (t) 


"1 = (1 + z ) 


1 

Wz) 


dH(z) 

dz 


- 1 . 


(3.19) 


With the use of Eq. (3.13) it follows that the deceleration parameter can be expressed as 


3 . . 2 k 2 + 3ct 4 — 4 a 5/7 uiC 2 

2 2 k 2 + 3 o 4 — 605/7 g 3 /7 2 ( 2 k 2 + 3o 4 — 605/7) 


(3.20) 


or, equivalently 
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(3.21) 


_ 3 1 + m- (2ni/3)H _w£i_ 

^ 2 1 + m — n\H a 3 H 2 (1 + m — n\H) 


As a function of redshift the deceleration parameter is obtained as 

3 1 + m — (2/3) nh(z) ws(l + z) 3 

2 U 1 + m — nh(z) h 2 (z) [1 + m — nh(z)] 


(3.22) 


The sign of the deceleration parameter indicates the nature of the expansionary evolution. 
If <7 > 0, the cosmological expansion is decelerating, while negative values of q indicate an 
accelerating dynamics. 

In the following we consider the cosmological implications of the present model for 
several equations of state of the cosmological matter. We first investigate the high density 
phase of the evolution of the universe with matter obeying the stiff and radiation equations 
of state. We then analyze in detail the behavior of dust (zero thermodynamic pressure) 
cosmological models. 


3.3 The stiff and radiation fluid cases 

We begin our study of the cosmology of the Lorentz violating Galileon model by analyzing the 
high density universe described by the Zeldovich (stiff) equation of state and by the radiation 
equation of state, respectively. The Zeldovich equation of state p = p is valid for densities 
significantly higher than nuclear densities p n , p > 10 p n . It can be obtained by constructing 
a relativistic Lagrangian which allows bare nucleons to interact attractively through scalar 
meson exchange, and repulsively through the exchange of a slightly more massive vector 
meson [21]. In the non-relativistic limit both the quantum and classical theories lead to 
Yukawa-type self-interaction potentials. But at the highest possible matter densities the 
vector meson exchange dominates. By using a mean field approximation for the nuclear 
interactions, it follows that in the extreme limit of very high densities the thermodynamic 
pressure tends to the energy density, p p. In this high density limit the speed of sound c s 
tends to one, (? s = dp/dp —>• 1. Therefore the Zeldovich (stiff fluid) equation of state satisfies 
the causality condition with the speed of sound less than or equal to the speed of light. The 
radiation fluid satisfies an equation of state of the form p = (1/3 )p. 

The variations of the Hubble parameter, energy density and deceleration parameter are 
represented, as functions of the redshift z, and for different values of the parameters m, n and 
s for a stiff fluid filled universe, in Figs. 1 and 2, while the variation of the same quantities 
for a radiation filled universe are represented in Figs. 3 and 4, respectively. 

In both the stiff and radiation fluid cases we have studied the evolution of the high 
density universe in the redshift range 5 < z < 20, with the initial condition h( 5) = 10. The 
behavior of the cosmological models for both equation of states is similar. For the range 
of considered parameters m, n and s, the Hubble function is a monotonically increasing 
function of 2 (monotonically decreasing in time), as is the energy density. The evolution 
of the deceleration parameter shows a decelerating expansion, with values of the order of 
q( 5) ~ 2 in the case of the stiff fluid universe and with q( 5) ~ 1 in the case of the radiation 
fluid. Of course, modifying the numerical values of m, n and s may lead to significantly 
different values of the cosmological parameters. 
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Figure 1. Variation of the Hubble parameter h(z) (left figure) and of the energy density r(z) (right 
figure) as a function of z for a stiff fluid filled universe, for different values of the parameters m , n , s: 
m = 0.0001, n = 0.0002, s = 0.0003 (solid curve), m = 0.0004, n = 0.0006, s = 0.0008 (dotted curve), 
to = 0.0006, n = 0.0008, s = 0.001 (short dashed curve), m = 0.0008, n = 0.001, s = 0.0012 (dashed 
curve), and to = 0.001, n = 0.0012, s = 0.0014, respectively. 



Figure 2. Variation with respect to z of the deceleration parameter q(z) for a stiff fluid filled universe, 
for different values of the parameters m, n, s: m = 0.0001, n = 0.0002, s = 0.0003 (solid curve), 
to = 0.0004, n = 0.0006, s = 0.0008 (dotted curve), m = 0.0006, n = 0.0008, s = 0.001 (short dashed 
curve), to = 0.0008, n = 0.001, s = 0.0012 (dashed curve), and to = 0.001, n = 0.0012, s = 0.0014, 
respectively. 

3.4 The universe filled with dust 

For a universe filled with dust, uj = 0, the time evolution equation Eq. (3.13) for the Hubble 
parameter takes the form 


(1 + m - mH)H + - 


1 + 777,-77 1 H 

3 


H 2 = 0. 


(3.23) 


Eq. (3.23) admits a de Sitter type solution H = Hq = constant, corresponding to 

3(1 + 777) 


Hn = 


2 ni 


(3.24) 


For this choice of parameters the expansion is exponential with the scale factor given 
by a = aoe Hot , with the deceleration parameter having the value q = —1. During a de Sitter 
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Figure 3. Variation as a function of z of the Hubble parameter h(z) (left figure) and of the mat¬ 
ter energy density r(z) (right figure) for a radiation fluid filled universe, for different values of the 
parameters m, n, s: m = 0.0001, n = 0.0002, s = 0.0003 (solid curve), m = 0.0004, n = 0.0006, 
s = 0.0008 (dotted curve), m = 0.0006, n = 0.0008, s = 0.01 (short dashed curve), m = 0.0008, 
n = 0.01, s = 0.012 (dashed curve), and m = 0.001, n = 0.012, s = 0.014, respectively. 



z 


Figure 4. Variation with respect to z of the deceleration parameter q(z) for a radiation fluid filled 
universe, for different values of the parameters m, n, s: m = 0.0001, n = 0.0002, s = 0.0003 (solid 
curve), m = 0.0004, n = 0.0006, s = 0.0008 (dotted curve), m = 0.0006, n = 0.0008, s = 0.01 (short 
dashed curve), m = 0.0008, n = 0.01, s = 0.012 (dashed curve), and m = 0.001, n = 0.012, s = 0.014, 
respectively. 


type phase the energy density of the universe is given by p = C 2 /a 3 , and it tends exponentially 
to zero. For arbitrary values of the parameters the general solution of Eq. (3.23) is given by 


3 , ,1 ni 

X (t — t 0 ) = — + ——-r 

2 H 3(1 + m) 


In ■ 


H 


3 + 2m — 2n\H 


(3.25) 


The variations with respect to the redshift 2 of the Hubble function, matter energy 
density and deceleration parameter for the dust universe are represented in Figs. 5 and 6 . 

Both the Hubble function and the matter energy density, shown in Figs. 5, are mono- 
tonically decreasing functions of time. The deceleration parameter, represented in Fig. 6 
has negative values in the range 0 < z < 2 , showing that for the given numerical values of 
parameters m, n and s, the universe experiences an accelerated expansion. 
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Figure 5. Variation as a function of z of the Hubble parameter h(z) (left figure) and of the matter 
energy density r(z) (right figure) for a dust universe, for m = 0 . 002 , s = 0 . 20 , and for different values 
of the parameters n: n = 1.653 (solid curve), n = 1.663 (dotted curve), n = 1.673 (short dashed 
curve), n = 1.683 (dashed curve), and n = 1.693, respectively. 



Figure 6. Variation as a function of z of the deceleration parameter q(z) for a dust universe, for 
m = 0.002, s = 0.20, and for different values of the parameters n: n = 1.653 (solid curve), n = 1.663 
(dotted curve), n = 1.673 (short dashed curve), n = 1.683 (dashed curve), and n = 1.693, respectively. 

4 The special matter-aether coupling cj)T 

In this Section we consider a special choice for the coupling between matter and the aether 
field. This coupling is common in massive gravity theories and galileons theories. In order 
to do this, we add a ae(j)T term to the action Eq. (2.1), where a q is a coupling constant, 
and T is the trace of the matter energy-momentum tensor. One can easily see that this new 
coupling will add a new term a^T to the left hand side of the aether equation of motion, 
Eq. (2.11), and a term of the form to the right hand side of the metric held equation 
Eq. (2.7), where 


V = «6 Ut iw + \(T- 2C m )d>g ia/ + , (4.1) 

with C m as the matter Lagrangian (For more detailed calculations, see [27]). 

We note that the energy-momentum tensor of the ordinary matter is not conserved 
due to the non-minimal interaction between scalar held and the energy-momentum tensor. 
In this case, the non-conservation equation can be obtained easily by taking the covariant 
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divergence of equation (2.7) and using equation (2.11), with the result 


V^T —_—_ 

V 1 + a 6 0 


~ + 0Vjy(£ m — —T) — 2 SI^[(j)B 


(4.2) 


where we have defined B fJU = g a P ■ The non-conservation of the energy-momentum 

tensor implies that the point particle does not follow the geodesic equation. In order to obtain 
the equation of motion for a point particle, we take assume that the energy-momentum tensor 
is described by a pressure-less perfect fluid, with 


I'nv — pUfj,u u . 

Substituting this into (4.2) and using the relation 


h vX A7 p T^ = u p A7 tx u x = ^ 


d 2 x x \ 

+ r A 


dx p dx a 
pu ds ds 


we obtain the equation of motion as 

d x P dx° 


cTar \ 

+ r A 


h vX 


Oi 6 


ds 2 ' pa ds ds p 1 + a§(j) 

In the case = 0 we obtain the standard geodesic equation 


-<j)V u p + pcf) p 


(4.3) 


(4.4) 


(4.5) 


4.1 The cosmological evolution in the presence of the 4>T coupling 

In order to obtain the effect of this new coupling term on the cosmological evolution of the 
universe, we note that one can choose the Lagrangian for a perfect fluid as C m = —p. Then 
for a flat Friedmann-Robertson-Walker geometry we obtain the metric and aether equations 
as 

— (15ct4 + 2 k 2 )H 2 — 21«5 H 3 — 3asH — p + A + -^(3 p — p)t = 0, (4.6) 

(3«4 + 2k 2 — Qa^H)H + -(3«4 + 2 k 2 )H 2 — 6 a$H 3 + p + -y(5p + p)t = 0, (4-7) 

and 

Q{lha^H 2 — l2a 4 H + a^)H + 90 a§H l — 108a 4 H 3 + I8013H 2 — 2{2>\H + A) — ae(3p — p) = 0, 

(4.8) 


where we assume that c\ = 0 in Eq. (3.3). 

As a first step in the analysis of the system of equations Eqs. (4.6) - (4.8) we rescale 
again the physical parameters according to 


«3 = —H 0 r], a 4 = — cr, a 5 = —H = H 0 h, 
o 15 ZIHq 

p = 2K 2 i7pr, p = 2 k 2 HqP, A = 2k 2 HqA, a§ = HqA, 


T 


t = 


Ho 


(4.9) 


where Hq is the present day value of the Hubble parameter and 77 , a and 6 are dimensionless 
constants. Then the system of Eqs. (4.6) - (4.8) takes the following dimensionless form 


|(1 + a)h 2 {r) - 9h 3 (r) — r]h(r) — r(r) + A(r) + y [3P(r) - r(r)] r = 0, (4.10) 
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dh(r) 

dr 



^0h 3 {r) + P + ^ [5 P{t) + r(r)] r = 0, 

(4.11) 


^0h 2 (r) - | <rh(r) + | 


*P- + ^V(r) - ™ah 3 (r) 

dr 7 v ' 5 w 


+ 6r//i 2 (r) — 2 


3A(t)/i(t) + 


dA(r) 

dr 


— A [3P(r) — r(r)] = 0. 


(4.12) 


The above equations can not be solved analytically. Hence in the following we consider their 
numerical solutions for two particular choices of the equation of state of the cosmological 
matter. 


4.2 The radiation fluid cosmological model 

In the case of a universe filled with a radiation fluid, P(t) = r(r)/3, the system of Eqs. (4.10) 
- (4.12) takes the form 

r(r) = |(1 + v)h 2 {r) - 0h 3 (r) - rjh{r ) + A(r), (4.13) 


dh (r) 1 

dr 21 (a + 5) — 309h(r) 


35(4At + l)A(r) + 35 /i(t)(4A?/t + r/) 


+ 5#(28Ar + 13)/i 3 (r) - 42 h 2 (r) [5A(<r + 1 )r + 2a + 5] , 


(4.14) 


IT = 105[7( g + 5)-100Mr)] { 35(4AT + 1)A(t) b 35 f - 75W>2 M + 

— 35/i(r)(4A?/T + rj) — 5#(28Ar + 13)/i 3 (t) + 42 h 2 (r) [5A(<r + 1 )r + 2a + 5] ^ 

ic: io 

+ 3 rjh 2 (r) + — #/i 4 (r) — 3/i(r)A(r)- —<t/i 3 (t). (4.15) 

( o 

The deceleration parameter can be obtained as 


q(r) = 


3h 2 (r)(lO0/i(r) — 7 (a + 5)) 


7 h(r) ^5(4A^r + 77 ) + 5/i 2 (t)(4A#t + 6) 


— 3/i(t)(10A(ct + 1 )r + 3cr + 5) I — 35(4Ar + l)A(r) 


(4.16) 


The time variations of the Hubble parameter, of the energy density, of the Lagrange multi¬ 
plier, and of the deceleration parameter of the radiation fluid filled universe are presented, 
for different values of the parameter A, 

in Figs. 7 and 8. In all cases a = 0.001, 9 = —0.001, and rj = 0.001. The initial 
conditions used to numerically integrate Eqs. (4.14) and (4.15) are h( 0) = 10 and A(0) = 0.5. 
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Figure 7. Variation with respect to r of the Hubble parameter h(r) (left figure) and of the energy 
density r(r) (right figure) for a radiation fluid filled universe, for different values of the parameters 
A: A = 1 (solid curve), A = 3, (dotted curve), A = 5 (short dashed curve), A = 7 (dashed curve), 


and A = 9 (long dashed curve), respectively. 
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Figure 8. Variation with respect to r of the Lagrange multiplier A (r) (left figure) and of the 
deceleration parameter q(r) (right figure) for a radiation fluid filled universe, for different values of 
the parameters A: A = 1 (solid curve), A = 3, (dotted curve), A = 5 (short dashed curve), A = 7 
(dashed curve), and A = 9 (long dashed curve), respectively. 


In the plots we have adopted constant values for the parameters a, 6 and 77 , and we have 
varied the numerical value of A, describing the strength of the coupling between cf) and T. 
As one can see from Fig. (7), the Hubble parameter is a monotonically decreasing function 
of time. The Lagrange multiplier A, presented in Fig. ( 8 ), has a complex cosmological 
evolution, consisting of two phases: a short, almost linearly decreasing with time which ends 
with reaching a minimum value at a finite time followed by a monotonic increase phase. The 
energy density r, shown in Fig. (7), is a monotonically decreasing function of time. As seen 
from Fig. 2.12, the dynamics of the deceleration parameter indicates that for all the chosen 
numerical values of the model parameters the radiation filled universe is in a decelerating 
phase. The time spent in the decelerating phase strongly depends on the numerical values 
of A. 

4.3 The dust cosmological model 

As a second example of a cosmological model in the framework of the SEA theory with 
aether-matter coupling we consider the case of a universe filled with dust, that is, matter 
with zero pressure. By taking P = 0 the cosmological equations describing the evolutionary 
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dynamics of the universe become 


= -r?/t(r) - dh 3 (r ) + (3/2)(<r + 1 )h 2 {r) + A(r) . . 

1 ’ Ar/2 + 1 ’ 1 ' ’ 

dh(r) —35ArA(r) + 35Ar/r/i(r) + 50(9Ar + 4)/i 3 (r) — 21L 2 (t) [A(3<r + 5 )r + a + 5] 
dr (At + 2) [7 (a + 5) — 10 0h{r)] 

(4.18) 


dA (r) 
dr 


1 


14AA(t)(— 5t]t + a + 5) 


2(Ar + 2) [7(a + 5) - 1O0 /i(t)] L 
+ /i 2 (r)(A(4ry(50 — 63crr) + 21 (<j + 1)(<t + 5)) + 300(4 — 3At)A(t) 

+ 42 rj(a + 5)) + 2h(T)(7Arj(5r]T — a — 5) — A(r)(A(lO0 + 21(5 — 3<r)r) + 42 (<r + 5))) 
+ 4/i 3 (t)(9At(5t?0 + 7 a (a + 1)) — A0(llcr + 25) — 20^0) + 15OA0 2 t/i 5 (t) 


+ 20/i 4 (t)( 1OA0 - 6A(32cr + 25)r + 9a - 75) 


(4.19) 


The deceleration parameter can be obtained as 

7/i(t)(5At/t — /i(t)(2A(4u + 5 )r — 5A0r/i(r) + a + 5)) — 35AtA(t) 


q{r) = 


(At + 2)/i 2 (t)[1O0/i(t) — 7(a + 5)] 


(4.20) 


The time variation of the Hubble parameter h(r), of the energy density r(r), of the Lagrange 
multiplier A(r) and of the deceleration parameter q(r) are presented in Figs. 9 and 10. 




Figure 9. Variation with respect to the dimensionless time r of the Hubble parameter h{r) (left 
figure) and of the matter energy density r(r) (right figure), for a dust universe, for different values 
of the parameters A: A = —0.75 (solid curve), A = —0.65, (dotted curve), A = —0.55 (short dashed 
curve), A = —0.45 (dashed curve), and A = —0.35 (long dashed curve), respectively. 

Here, in a similar fashion to that of the radiation fluid filled universe model, we have also 
fixed the numerical values of the free parameters a, 0 and tj as a = 0.001, 0 = —0.005, tj = 
0.001, and have varied the coupling parameter A < 0. The initial conditions used to integrate 
the field equations Eqs (4.18) and (4.19) are h( 0) = 2 and A(0) = 1.5, respectively. The 
Hubble parameter, depicted in Fig. 9, is a monotonically decreasing function of time, as well 
as the energy density, shown in Fig. 9. The Lagrange multiplier A, plotted in Fig. 10, shows 
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Figure 10. Variation with respect to the dimensionless time r of the Lagrange multiplier A(r) (left 
figure) and of the deceleration parameter q(r) (right figure) for a dust universe, for different values 
of the parameters A: A = —0.75 (solid curve), A = —0.65, (dotted curve), A = —0.55 (short dashed 
curve), A = —0.45 (dashed curve), and A = —0.35 (long dashed curve), respectively. 


a similar behavior as for the radiation fluid cosmological model, decreasing monotonically 
with time. The deceleration parameter, depicted in Fig. 10, starts its evolution at a high 
redshift (high value of h) with a positive value of the order of q ~ 1, indicating an initially 
decelerating cosmological evolution. Due to the presence of the aether field and of its coupling 
with matter the deceleration parameter monotonically decreases in time and the dust filled 
universe enters an accelerating stage, with the deceleration parameter reaching values of the 
order of q ~ — 1. 

5 Discussion and final remarks 

In this paper we have studied a Lorentz violating theory of gravity by introducing a time- 
like vector field in the gravitational action. The time-like vector field is constructed from 
the convariant derivative of the scalar field <p. In order to impose the time-like property of 
the vector we have added the constraint + 1 = 0 to the action through a 

Lagrange multiplier. In order to add a kinetic term for the scalar field <f>, we note that the 
canonical kinetic term for the scalar field is already used in the constraint equation. With 
this constraint term, the scalar field becomes dynamical in the theory [26] . In order to enrich 
the dynamics of the scalar field one may add some higher order derivative scalar interaction 
terms to the action. This was done in [13] where the authors used the most general kinetic 
terms of a vector field to construct such interaction terms. However, the resulting theory 
suffers from ghost instabilities because of these higher order derivative interactions. This 
can be seen by noting that the theory [13] is equivalent to the projectable version of Horava- 
Lifshitz gravity which suffers from instabilities and strong coupling in low energies [7, 13, 16]. 
In this paper we have added some special higher order derivative interactions which produce 
second order field equation, and hence, no Ostrogradski instabilities. These interactions are 
known as Galileon terms [17]. 

The present model can also be considered as a generalization of the mimetic dark mat¬ 
ter models [23]. In the mimetic dark matter model, one defines a physical metric out of a 
scalar field and some primary metric tensor in such a way that the physical metric becomes 
conformally invariant. Thus constructing a gravitational theory by the physical metric will 
automatically produce a conformally invariant gravitational theory. However, as was men- 
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tioned in [24], this theory is equivalent to a gravitational theory coupled to a scalar field 
and subject to the constraint that is tinre-like. This is in fact equivalent to the scalar 
Einstein-aether theory [13]. In this sense, our present model is a generalization of the mimetic 
dark matter theory with some additional healthy derivative self-interaction terms. This will 
actually make the dark matter imperfect [25]. The problem of whether this model can satisfy 
dark matter experimental data will be considered elsewhere. 

In the present paper, as a first step in the in-depth investigation of the model, we 
have studied the cosmological consequences of the Lorentz violating Galileon theory. If one 
assumes that the matter content of the universe consists of radiation or stiff matter, then the 
theory predicts a decelerating universe. In this case the universe decelerates more rapidly 
with stiff matter as compared to radiation. For a dust baryonic matter the model suggests 
an accelerating universe with an exponential de Sitter acceleration at late times. 

We have also considered a special coupling between ordinary matter and the scalar 
field. This coupling will break the conservation law of the ordinary matter and hence one 
can expect some modification in the predictions of the Solar System tests, which may put a 
constraint on the value of the coupling constant ag in Eq. (4.1). This will also be done in a 
separate work. In this paper, we have considered only the cosmological implications of this 
extra term. In the radiation or stiff matter cases, the quantitative behavior of the universe 
does not change as compared to previous models without the interaction term, predicting 
a decelerating universe. In the case of a dust baryonic matter, however, one can obtain 
a universe starting from a decelerating phase and entering into an accelerating phase at 
later times. Other astrophysical and theoretical implications of this model such as black hole 
solutions and the stability of the cosmological solutions will be studied in a future publication. 
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